2D compressibility of surface states on 3D topological insulators 
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We develop a theory for the compressibility of the surface states of 3D topological insulators and 
propose that surface probes of the compressibility via scanning single electron transistor microscopy 
will be a straightforward way to access the topological states without interference from the bulk 
states. We describe the single-particle nature of the surface states taking into account an accurate 
Hamiltonian for the bands and then include the contribution from electron-electron interactions 
and discuss the implications of the ultra-violet cutoff, including the universality of the exchange 
contribution when expressed in dimensionless units. We also compare the theory with experimentally 
obtained ^ as extracted from angle-resolved photoemission spectroscopy measurements. Finally, 
we point out that interaction-driven renormalization of the Fermi velocity may be discernible via 
this technique. 



Recently, there has been great interest in trying to en- 
gineer materials with topologically protected states be- 
cause of the apphcations that such states may find in 
quantum computing. In particular, it was predicted that 
certain 3D layered crystals with strong spin-orbit cou- 
pling and an inverted conduction band would have an 
insulating band gap in the bulk, but topologically ro- 
bust conducting surface states with linear low-energy 
dispersion protected by time- reversal symmetry.^ Angle- 
resolved photo-emission spectroscopy (ARPES) experi- 
ments on crystals such as Bi2Te3 and Bi2Se3 have proven 
the existence of these mid-gap surface states with nearly 
linear spectra. However, in most current samples the 
Fermi energy is well above the bulk gap meaning that 
transport measurements of the surface states are difficult 
because current can also flow in the bulk of the crystal®, 
i.e. the bulk is a conductor rather than an insulator even 
at low temperature^, although this can be mitigated to 
some extent by using thin films of the topological insula- 
tor (TI) crystal® or by doping the crystal to compensate 
for the bulk conduction)^ Therefore, it is essential that 
methods are found to both control and characterize these 
topologically protected states. 

In this Rapid Communication, we propose using sin- 
gle electron transistor (SET) microscop y^"' as a way of 
mitigating the effect of the bulk states by approaching the 
surface directly. Local capacitance measurements made 
via this technique can be converted straightforwardly into 
the quantity ^ (where /i is the chemical potential and n 
is the carrier density in the surface state) which is directly 
linked to the band structure and to the electronic com- 
pressibility. A great advantage of measuring this ther- 
modynamic compressibility is that it incorporates single- 
particle many-body renormalization directly,^^ and thus, 
quantitative information both about the single-particle 
(SP) band structure and many-body renormalization are 
obtained in one stroke as a function of 2D carrier den- 
sity (or equivalently, the 2D Fermi energy of the surface 
states). We begin by calculating ^ for two different 
commonly used approximations for the noninteracting 
band structure of the Hamiltonian for the surface states 
of TIs such as Bi2Se3 and Bi2Te3 and then take into 



account the electron-electron interaction demonstrating 
that this contribution is universal and identical for both 
band structure approximations when expressed in dimen- 
sionless units. However, we demonstrate that the high 
dielectric constant of these materials reduces the quan- 
titative effect of the interactions. We also argue that it 
may be possible to observe renormalization of the quasi- 
particle Fermi velocity at sufficiently low density by this 
technique. We then compare the results of our analyti- 
cal SP calculations to estimations of ^ extracted from 
ARPES measurements of the band structure. Finally, we 
comment on the role of disorder and charge inhomogene- 
ity on ^ in these systems. 

If the plane of the surface corresponds to the (111) 
crystal direction {i.e. parallel to the layered structure 
of the lattice) then the simplest approximation for the 
band structure is the commonly-used linear approxima- 
tion given by the Hamiltonian^ ij' = hvpa ■ k where 
is the band velocity, k is the two-dimensional wave 
vector in the plane of the surface, and <yx,y sue Pauli 
matrices in the real spin space. The spectrum of this 
Hamiltonian is trivially the linear Dirac dispersion given 
by e^j, = Xhvpk where A = ±1 denotes the band. The 
corresponding eigenvectors are two-component spinors in 
the spin space. The 2D density of carriers in the SP limit 
is given by n = kp/{ATr) where kp is the Fermi wave vec- 
tor and as a matter of convention, we say that zero den- 
sity corresponds to the chemical potential being located 
at the Dirac point. Elementary manipulations show that 
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A more sophisticated approximation for the band 
Hamiltonian including the curvature in the valence band 
which is noticeable in the ARPES data is H'^ = ^-4r + 
fivpa ■ k where m* is the effective mass of the carriers 
The spectrum associated with this Hamiltonian is e\^, = 

+ Xhvpk. These spectra are shown in Fig. [TJa). 
In the case with the quadratic correction, there is an 
unphysical turning point in the valence band which de- 
mands that some care must be taken when describing this 
band. The turning point is located at at fco = vpm* /h, 
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FIG. 1. (a) Single particle bands. The dotted line is the 
'turned' band structure where the quadratic band is inflected 
at its turning point, (b) Single particle ^ corresponding to 
Eqs. m and (111). 
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FIG. 2. (a) The self-energy relative to its value at zero density. 
Solid lines indicate kc = 1/a and dashed lines show kc — 2-k ja. 
(b) The self-energy normalized by its absolute value at zero 
density. Lines for both dispersions and all values of fee collapse 
onto each other. Throughout, the linear band is shown by red 
lines and the quadratic band by black lines. 



energy contribution from the exchange interaction gives 

Scx,A(k) = -^ y d2q^T/c(q)/A'(k-Hq)x 

X [1 + AA' cos(0fe+g - dk)\ (3) 

where VcCq) = 27re^/(K|q|) is the two-dimensional 
Fourier-transformed unscreened Coulomb interaction, k 
is the static dielectric constant of the environment, /A(k) 
is the occupancy of the state with wave vector k in band 
A, 9k is the angle that the wave vector k makes with 
the X axis. It is well-known that for the compressibility 
the exchange energy (or equivalently, the Hartree-Fock 
term) is the dominant interaction correction, and this 
is even more true for TI systems because of their very 
large background dielectric constant. Therefore we antic- 
ipate very small interaction corrections beyond exchange 
in this problem. Note that the only place that the band 
structure enters into this expression is in the occupancy 
factors so that it is applicable to systems described by 
both and H'^ . At zero temperature the angular part 
of the integration of q can be computed analytically and 
yields 

where the appropriate SP energy must be inserted into 
the step function. 
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!^K{X)^{k' + k)E{X) AA' = -1 



the energy at this wave vector is Sq 



and the as- 



sociated density is no 
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/{Anh ). Hence, defin- 



ing wave vectors, energies, or densities which are greater 
in magnitude than fco, eo, and no requires a specific def- 
inition of the band structure, and this will become vital 
when the interaction effects are incorporated. We defer 
further discussion of this specific point until it is rele- 
vant. However, for n > no we use the relationship for 
the density above, and find 



dfi 2nsgn{n)h 
dn 



m 



(2) 



where sgn(a;) denotes the sign of the argument. Figure 
[ijb) shows the single-particle ^ for each approximation 
as a function of density. Unless otherwise stated, we use 
vf = SxlO^ms^^, m* — Q.2me {nie is the electron mass), 
and a — 0.41nm throughout this paper. The quadratic 
correction preserves the dominant 1/Vn behavior at low 
density, but introduces an asymmetry between the con- 
duction and valence bands. In addition, |^ goes to zero 
at no- 

We now include the effect of the electron-electron in- 
teraction on the compressibility. Evaluating the self- 



K and E are complete elliptic integrals of the first and 
second kind, and X = 2^/kk' /{k + k'). The radial in- 
tegral must be evaluated numerically. We note that the 
exchange integral of Eq. ([3]) has an ultraviolet high mo- 
mentum divergence arising from the linear dispersion, 
which must be regularized through a high wave vector 
cut off kc- As is usual in condensed matter physics, there 
is a real cut off in the momentum arising from the lattice 
structure, and therefore, kc ~ 1/a, and the interaction 
strength depends explicitly on the short-distance lattice 
cut-off in the theory. To illustrate the qualitative fea- 
tures of the physics, and to determine the dependence of 
^ on kc we start by describing the k = 1 case in which 
interactions are the strongest. This function is shown in 
Fig. [2ja) for both the linear band and the band with a 
quadratic correction for two different values of the ultra- 
violet cutoff. Details of the definition of the band struc- 
ture for the quadratic correction are given below. In Fig. 
[SJb), the same data is shown but the units are scaled to 
demonstrate how the self-energy depends on the cutoff 
kc- The linear and quadratic bands give identical results 
when the self-energy is scaled by its value at zero density. 
We emphasize that this result shows the universality of 
the exchange contribution, which is reminiscent of simi- 
lar results found in semiconductor heterostructures.— 



3 




FIG. 3. Linear band, (a) ^ as a function of density for the 
single-particle band, and two values of kc- This illustrates 
that the choice of kc is crucial: small kc gives a negative 
contribution to ^ at higher densities, large kc gives a positive 
contribution, (b) ^ as a function of log(fecffl) for different 
carrier densities. 

The independence of the dimensionless many-body cor- 
rections of the details of the TI band structure is an im- 
portant new result of our work. 

The exchange self-energy is incorporated into the cal- 
culation for ^ in the following way. The Hartree-Fock 
chemical potential is the sum of the SP kinetic energy and 
the exchange self-energy: /j,hf — ^J■ + '^cx and the cor- 
responding compressibility is given by d^np/dn, which 
may be computed numerically. When applied to the lin- 
ear band structure associated with H\ the expression 
in Eq. ([3]) contains an ultra-violet divergence and so the 
value of the high-frequency cutoff kc becomes important, 
as hinted at in Fig. [2] In Fig. [Sja) we show ^ for the 
linear band with the exchange contribution for two dif- 
ferent physically reasonable values of kc- In all cases, the 
n~2 behavior persists at low density. For kc = 1/a the 
exchange causes a reduction in ^ relative to the SP case 
for all but the lowest densities and may become negative 
for large valence band doping. In contrast, kc = 27r/a 
gives an enhancement to ^ . This indicates that the pre- 
cise value of the ultraviolet kc has an important quan- 
titative effect on In Fig. Mjo) we show the depen- 
dence of ^ on kc for three different values of the density. 
Noting the logarithmic scale on the kc axis, we see that 
^ oc log(fcca). Hence, in order to make quantitative 
predictions of the behavior of ^ at finite density, some 
physical intuition must be used for choosing the 'correct' 
value of kc- On the other hand, the logarithmic depen- 
dence of the interaction effect on kcU translates basically 
into a slow logarithmic renormalization of the coupling 
strength in units of kc/kp which could manifest itself at 
very low carrier densities.— 

For the quadratic band, Eq. ^ still holds, but the ap- 
propriate dispersion must be used in the Fermi function. 
The shape of the valence band introduces a further com- 
plication. Taken at face value, the band minimum at k^ 
should be interpreted as producing a Fermi surface with 
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FIG. 4. The quadratic correction, (a) ^ for kc^ko. (b) ^ 
for various values of kc > ko. 

non-trivial topology, and contributions to ^ should be 
included from both branches of the valence band. How- 
ever, this does not correspond to the experimental data 
(see, for example, Ref. |l7l ) where the turning point is 
not seen. The first (and most intuitive choice) is to set 
kc = ko- This may be physically reasonable because the 
topological surface bands typically merge with the bulk 
bands at roughly this wave vector. Another possibility 
is to artificially turn the band structure for \k\ > kg in 
the valence band (as illustrated by the dotted line in Fig. 

[T^a)) so that for |fc| > ko we have e^^, = — |;;|r — hvpk. 
Then, from a technical point of view, kc can be set ar- 
bitrarily high and the quadratically-dispersing nature of 
the band is retained. Figure IH^a) shows the first of these 
cases. Since kc is relatively small, the density can reach 
the regime where the self-energy is decreasing rapidly, 
n/n{kc) ~ 0.1. In that case, the negative slope of the ex- 
change energy is stronger than the positive slope of the 
SP chemical potential so that /iHFA is a decreasing func- 
tion of n for n < implying that cI^yifk/ dn is negative. 
We believe that some of these interesting questions about 
the precise value of the ultraviolet cut-off in this problem 
can actually be answered through a careful comparison 
between experiment and theory on the TI compressibil- 
ity. 

For the turned band structure with kc > fco, the evolu- 
tion of Sex with increasing kc means that the sign of the 
exchange contribution can be positive or negative. For a 
given density, we have |^ oc log(A:ca) as was the case for 
the linear band. 

The static dielectric constant of crystals such as Bi2Se3 
and Bi2Te3 is estimated to greater than 50 which 

will naturally have the effect of reducing the strength of 
the interactions. If we assume that the effective dielectric 
constant is the average of that in the material and that of 
the air then we can take a physically reasonable value of 
K = 20. Figure O shows that the exchange contribution to 
1^ is reduced so much that it is essentially non-existent. 
Hence, for the purposes of examining the compressibility 
of the surface states, the SP calculation is likely to be 
quantitatively sufficient. 

At this point, we pause to comment on another ex- 
perimental manifestation of the exchange self-energy. In 
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principle, the Fermi velocity is renormalized by this 
interactioniSi^ and there has been some experimental 
hint of this in graphene in both zero magnetic field^i and 
high magnetic fieldsi^"— The renormalized quasiparticle 
velocity can be written as 



1 dfi 



HF 



kp 



with Ts = I (JiKVp) and thus the interaction-induced 
correction will be proportional to log(A:c/fci?). Of course, 
the Fermi velocity measured in experiment includes this 
renormalization so it is not directly measurable, but since 
the self-energy is a non-monotonic function of the wave 
vector, its effect may be seen by a density-dependent de- 
viation of the Fermi velocity from a constant value which 
represents the combination of the SP and mean interac- 
tion contributions. Such a nonlinear reconstruction of 
the 2D TI Fermi surface at low carrier density will be a 
direct manifestation of its Dirac spectrum and the asso- 
ciated ultraviolet renormalization effect familiar in quan- 
tum electrodynamics. 

We can also extract ^ from ARPES measurements of 

an 

the band structure. Digitizing curves from Zhu et alh 
gives the band structure of Bi2Se3 in the KVK direction 
from which we can numerically extract ^ for compari- 
son to theory. However, this process is complicated by 
the fact that the curves on either side of the F point are 
not perfectly symmetric, it is also known that there is a 
hexagonal distortion to the Fermi surface^? and ARPES 
data contains some noise. To combat these difficulties we 
apply a Gaussian smoothing convolution to the experi- 
mental data with a width a = O.OSnm^^ to both branches 
shown in the KTK plot and take the average of the two 
to arrive at the approximate band structure. Then we 
associate the Fermi wave vector kp with this averaged 
curve for each value of energy and from this compute the 
density as n = kp/{A'!T) which assumes a circular Fermi 
surface. A numerical derivative of the chemical poten- 
tial with respect to the density may then be taken. The 
results of this procedure are shown in Fig. Eljb) along 
with a least-squares best fit of the quadratic correction 
dispersion using m* and vp as fitting parameters from 
the conduction band data. This fitting procedure yielded 
TO* = 0.417TOe and vp = 1.00 x 10''ms~\ both of which 
are larger than the current common estimates for these 
parameters in Bi2Se3. The dashed line corresponds to 
the predicted ^ for a purely quadratic band, which is 

constant at a value of ^ = « l.lSeVnm^. This in- 

an m 

dicates that the quadratic part of the band is completely 
dominated by the linearity of the band structure over 
the whole experimentally pertinent range of carrier den- 
sity very clearly establishing that the Dirac spectrum is 
dominating the physics in these systems. 

One likely challenge to using the compressibility to in- 
vestigate the topological surface states is the observed 
existence of substantial disorder-induced inhomogeneity 
in the charge landscape of the surface.^ Surface probes 
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FIG. 5. (a) The effect of electron-electron interactions for the 
quadratic band with «: = 20 and kc = 1/a. (b) Black line: ^ 
computed from band structure extracted from ARPES data 
in Zhu et ali^ (see text for details). Red line: least-squares 
best fit to the analytical result of Eq. ((2]). The parameters 
extracted from the conduction band are m* = 0.417me and 
VF = 1.00 X 10^ms"\ 



such as scanning tunneling microscopy have shown that 
the screening of an external disorder potential created 
by charged impurities in the lattice of Bi2Te3 and Bi2Se3 
lead to the formation of 'puddles' of electron and holes 
with spatial extent of the order of lOnm and associated 
fluctuations in the local chemical potential of approxi- 
mately lOmeV. This fluctuation has previously been in- 
vestigated in the context of charge transport?^ but in 
the case of our proposed SET measurements, it should 
be pointed out that the area of the sample which influ- 
ences the tip has a radiusiSiii of the order of lOOnm, 
much larger than one puddle. Therefore, the SET will 
experience some averaged field corresponding to the in- 
homogeneous landscape. The effect of the inhomogeneity 
is that even when the global average of the density is zero, 
the local density is always finite and may be positive or 
negative. Therefore, in the inhomogenous case, the lo- 
cal density never approaches zero and hence the ^/^/n 
divergence is curtailed. At high density, the slow change 
of ^ with density implies that the inhomogeneity will 
have little effect on the measured compressibility. There- 
fore, the extraction of band parameters from SET data 
is likely to be most accurate at higher density. 

In this Rapid Communication, we have proposed that 
a capacitive surface probe such as SET microscopy will 
reveal the properties of the topological surface states of 
crystals such as Bi2Te3 and Bi2Se3 via measurement of 
the compressibility, or 4^. We have analytically com- 



puted the SP contribution to ^ for two different model 
band structures and described the role that electron- 
electron interactions and charge inhomogeneity in the 
crystal will play. In summary, the high dielectric con- 
stant implies that the quantitative effect of interactions 
is negligable (although at low density there is some hope 
of detecting a renormalization of the Fermi velocity due 
to electron-electron interactions) , and the inhomogeneity 
can be accounted for via a straightforward phenomeno- 
logical averaging technique which will result in the low- 
density divergence being cut off. We have established 
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that compressibility measurements, when compared with SP and many-body properties of 2D TI surface states, 
theory, could provide valuable information about both We thank US-ONR for support. 
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